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Department of Physics, Zhengzhou University, China. 

1 . Introduction 

The earth's atmosphere is continually being bombarded by primary 

cosmic ray particles which are generally believed to be high-energy 

nuclei. The fact that the majority of cosmic ray primaries are charged 

particles and that space is permeated with random magnetic fields, 

means that the particles do not travel in straight lines. This makes 

the identification of distant sources very difficult. Nevertheless, 

studies of the arrival time and direction distribution of cosmic rays 

are still used to seek significant information on problems of their 

origin. Prom the beginning of the 1950's to the middle of the 1960's 

about 50 experiments were carried out to study the arrival directions 

tiu 17 

of EAS with energy range from about 10 to 10 eV. Lins ley and 
Watson (1 ) summarised the results of these experiments and gave a 
review at 15th Cosmic Ray Conferees. 

On the other hand, the arrival time distribution of EAS may 
also transfer some information about the primary particles. Actually, 
if the particles come to our, earth in a completely random process, 
the arrival time distribution of pairs of successive particles should 
fit an exponential law. This is derivable from Poisson's distribution. 
Several groups £ 2,5 5 have reported a non-random component in the 
arrival times of EAS with E>10 ,<f eV, but others(4} <iid not find it. 

The work reported here was carried out at Sydney University 
from May 1982 to January 1983. The results are discussed and compared 
with that of some other experiments. 

2. Experiment System 

This experiment work was carried out by using the Sydney 
Small Air Shower Array. This sea level array has been described in 
detail elsewhere^5j. The array is composed of four fast timing 
scintillators arranged as a 25 mjf 25 m horizontal square with four 
triggering scintillators in a 4m)f4m square in the same plane as the 
fast timing scintillators. The triggering scintillator square is close 
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but not at the centre of the fast timing square. Pulses from the 
scintillators are sent via fast rise time coaxial cables to the 
recording electronics which is in an air conditioned hut with the 
triggering scintillators. The recording electronics uses NIM/CAMAC 
System. A master trigger pulse is produced if three out of four of 
the triggering scintillators receive shower particles. Showers are 
analysed in a desk-top compuer which is employed to control the 
recording system. The shower direction is found by a weighted least- 
squares fit of a plane shower front to relative arrival time of each 
detector. Then, using the maximum likelihood method to estimate the 
shower size and core position. Finally, the celestial coordinates of 
the incident direction of the event is found out. The rate of analysed 
showers is 127 per day or one every 11.5 minutes on average. The tine 
needed for the full analysis is about 50 seconds. 

5. Analysis and results 

The data discussed here were obtained in the period from May 27 
in 1982 to January 12 in 1985 at the campus of Sydney University 
(latitude 53 *)• During this period more than 25000 showers were 
recorded. The data are analysed in two different ways as follows. 



Fig. 1 Cosmic ray arrival time distribution 


Firstly, all the data are used to plot the arrival— time 
distribution of the events, that is, the distribution of time-separation 
t between consecutive events on a 1 minute bin size (Fig. 1 ). The smooth 
curve shows the expected exponential distribution of the arrival times 
assuming that whs time of occurrence of the events is completely randoms 

n-N«(e“ Xt '- e” Ktl ) 

where x is the average number of events per unit time. As can be seen 
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from Pig. 1, the observation data are compatible with random 
expectation. After adding the estimated losses for the observation 
data, no deviations are greater than 3<T. So that, no experimental 
evidence for abnormal behaviour in the inter shower arrival time 
distribution has been found. 

Then, the data are analysed with respect to the sidereal time 
variation. Since the experiment was interrupted occasionally for 
maintenance and by power failures, some allowances for this 
interruption must be made before the analysing. After rejecting 
certain numbers of events from those "over-exposure" time intervals, 
the run time for every sidereal hour interval is unified. The data 
are analysed by using the "random walk" harmonic method and the results 
are : 

the fractional amplitude r * %* 

the probability of observing an amplitude p = 0.21, 

the phase of maximum f * 74* ± 33° • 

To compare these results with that of the experiments summarised 
in Ref. (1}, they are plotted together in Pig. 2 (after Kiraly et al 
C63 ). the value of p of the present experiment sits at the middle of 
Pig. 2 (l), and the value of ir just falls in the gap of the old data 
in Pig. 2 (2). If anything this confirms the isotropy of arrival 
directions. 




Pig. 2 Comparison of the present experiment with the data of Ref. (1} . 
The ordinates are the frequencies and the circles contain 
the numbers of experimental sets contributing to the data 
(excluding the present experiment). 


The data are also tabulated on Mercator Projections and a series 
of chi-squared tests is applied. The results are shown in Fig. 3* The 
largest value of X 2 for the tabulation is 57 • The probability that X* 
should have exceeded this value for one or more of eight independent 



261 


OG 5.3 - 8 


Mercator projection of Celestial Sphere 
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Fig. 3 Cellular division of Mercator Projection of 
celestial sphere for 17164 showers. 


hands is 12%. Thus the chi-squared tests show no significant evidence 
of anisotropy in the tabulated distributions. 
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